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MULTI-SEGMENT CONE-BEAM RECONSTRUCTION SYSTEM AND METHOD 

FOR TOMOSYNTHESIS IMAGING 

5 RELATED APPLICATIONS 

This application claims priority to United States Provisional Application No. 
60/526,719, entitled Multi-Segment Cone-Beam Reconstruction Algorithm for 
Tomosynthesis Imaging and filed on December 3, 2004, which application is 
incorporated herein by reference.FIELD OF THE INVENTION 

10 The present invention relates to a system and method for imaging a target 

element using tomosynthesis. More specifically, the invention relates to a system, 
method and computer program product for creating a three-dimensional image of target 
elements from a plurality of radiation absorbance projection images taken from different 
angles. 

15 

BACKGROUND OF THE INVENTION 

Imaging of a patient's tissue has become a common screening and/or diagnostic 
tool in modern medicine. One example of such imaging is mammography, or the 
imaging of a patient's breast tissue. Breast cancer remains the most common cancer 

20 among women today, however, at this time there is no certain way to prevent breast 
cancer and the best strategy for dealing with breast cancer is early detection of the 
cancer so that it may be treated prior to metastatic spread. Accordingly, it is important 
for patients to have access to imaging techniques and systems that will detect very small 
cancers as early in their development as possible. 

25 A three-dimensional imaging approach called "tomosynthesis" has been 

developed (see United States Patent No. 5,872,828, which is incorporated herein by 
reference for its teachings relating to tomosynthesis systems and methods) which shows 
great promise for early detection of cancer. Tomosynthesis allows the reconstruction of 
a true volumetric distribution of absorption coefficients on the basis of the information 

30 contained in a series of projections acquired from a series of viewpoints about the target 
object. The viewpoint need not be regularly spaced, numerous, or arranged in any 
regular geometry. The tomosynthesis technique has been demonstrated to provide 
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useful spatial differentiation of overlapping and nearby tissues at \6ry high resolution 
comparable to projection 2D imaging, with approximately comparable radiation dose. 
The problem of 3D reconstruction from tomosynthesis projections is difficult 
5 one. One promising technique for 3D reconstruction from tomosynthesis projections is 
provided in United States Provisional Patent Application Serial No. 60/446,784, which 
is incorporated herein by reference. This technique applies a cone-beam geometry in an 
iterative forward-projection and back-projection method based on maximum-likelihood 
estimation of volumetric distribution of attenuation coefficients, using an estimation- 
10 maximization algorithm. However, the amount of computing power required to 

perform the 3D reconstruction will, for typical inexpensive computer systems, result in 
significant delays before the 3D reconstruction is available. 

SUMMARY OF THE INVENTION 

15 The present invention improves upon systems and methods known in the art by 

providing tomosynthesis apparatus and techniques for three-dimensional imaging of 
target elements that overcome the problems of conventional three-dimensional imaging 
systems and enable the use of tomosynthesis to efficiently provide accurate three- 
dimensional imaging of a target element in a shorter time than has previously been 

20 possible. The invention can increase the speed of the overall volume reconstruction by 
applying a segmentation of the target volume, while providing a quality of 
reconstruction that is substantially the same as for reconstructions that do not employ the 
segmentation algorithm of the invention. In preferred embodiments, the invention 
directs the reconstruction of volume segments to a multiplicity of CPUs within a multi- 

25 CPU computer cluster so that the volume segments can be reconstructed simultaneously 
in a parallel architecture. In further preferred embodiments, the volume segments are 
selected to be optimally overlapping in order to provide mutual boundary coverage of 
what would otherwise be under-determined edge volume segments. 

In one aspect, the invention provides a tomosynthesis method for creating a 

30 three-dimensional reconstruction of a target element volume through the acquisition of 
radiation absorbance projection images of the target element volume through a limited 
plurality of viewpoints. The target element volume is divided into a plurality of volume 
segments and a reconstruction algorithm is applied to each segment to generate a three- 
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dimensional reconstruction of each volume segment. The three-dimensional 
reconstruction of each volume segment is then merged to create a full three-dimensional 
reconstruction of the target volume. 
5 In a further aspect of the invention, a system for three-dimensional 

tomosynthesis imaging of a target element is provided having an image acquisition 
element and a processor. The image acquisition element obtains a plurality of images of 
the target element from a plurality of viewpoints. The image acquisition element 
includes a radiation source positionable at a plurality of viewpoints with respect to the 

10 target element and a radiation detector. The radiation detector is positioned so as to 

detect radiation emitted by the radiation source passing through the target element and 
determine a plurality of attenuation values for radiation passing through the target 
element to establish a radiation absorbance projection image of the target element for a 
particular radiation source position. The processor is configured to apply a 

15 reconstruction algorithm to the radiation absorbance projection images of the target 
element obtained from a plurality of radiation source positions to generate a three- 
dimensional reconstruction of the target element. The processor is further configured to 
divide the target volume into a plurality of image reconstruction volume segments for 
separate image reconstruction of the volume segments and to merge the reconstructed 

20 volume segments into a three-dimensional reconstruction of the target element. 

In a still further aspect, the invention provides a computer program for three- 
dimensional tomosynthesis imaging of a target element volume from a plurality of 
radiation absorbance projection images obtained at different angles from an image 
acquisition element having a radiation source positionable at a plurality of viewpoints 

25 with respect to the target element and a radiation detector positioned so as to detect 

radiation emitted by the radiation source passing through the target element volume and 
determine a plurality of attenuation value for radiation passing through the target 
element to establish a radiation absorbance projection image of the target element 
volume for a particular radiation source position. The computer program code is 

30 embodied in a computer readable medium and includes computer program code for 
dividing the target element volume into a plurality of volume segments, applying a 
reconstruction algorithm to each segment to generate a three-dimensional reconstruction 
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of each volume segment, and merging the three-dimensional reconstruction of each 
volume segment to create a three-dimensional reconstruction of the target volume. 

5 BRIEF DESCRIPTION OF THE DRAWINGS 

The invention will be more fully understood from the following detailed 
description taken in conjunction with the accompanying drawings, in which: 

Figure 1 A provides a schematic of a tomosynthesis system according to one 
10 aspect of the invention; 

Figures IB and 1C illustrate two orthogonal views of the geometry of the 
tomosynthesis system of Figure 1A with the view of Figure IB being along a patient's 

i 

chest wall (with the X-ray source traveling into and out of the page) and the view of J 

j 

15 Figure 1C being in a direction toward the patient's chest wall; 

Figure 2 illustrates segmented projections and volume coverage for the : 

i 

tomosynthesis system of Figure 1A and from the perspective of Figure IB; [ 

20 Figures 3A and 3B illustrate, respectively, non-overlapping and overlapping 

segmentation geometries that may be implemented in a method and system of the 
invention; 

Figure 4 illustrates diagrammatically a segmentation of a Tomosynthesis 
25 projection that may be implemented in a method and system of the invention; 

Figure 5 illustrates diagrammatically a reconstruction segmentation that may be 
implemented in a method and system of the invention; 

30 Figure 6 illustrates diagrammatically a reconstruction of a volume segmentation 

that may be implemented in a method and system of the invention; 
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Figure 7 illustrates diagrammatically a mergence of volume segmentations that 
may be implemented in a method and system of the invention; 

Figures 8 A, 8B and 8C illustrate a reconstructed volume image slice at Z=10 mm 
for a convention reconstruction, a segmented reconstruction of the invention, and a 
difference between the two, respectively; 

Figures 9A, 9B and 9C illustrate a reconstructed volume image slice at Z=30 mm 
for a convention reconstruction, a segmented reconstruction of the invention, and a 
difference between the two, respectively; 

Figures 10A, 10B and 10C illustrate a reconstructed volume image slice at Z=50 
mm for a convention reconstruction, a segmented reconstruction of the invention, and a 
difference between the two, respectively; and 

Figures 1 1 A, 1 IB and 1 1C illustrate a reconstructed volume image slice at Z=70 
mm for a convention reconstruction, a segmented reconstruction of the invention, and a 
difference between the two, respectively. 

DETAILED DESCRIPTION 

The systems and methods of the present invention improve upon systems and 
methods known in the art by providing tomosynthesis apparatus and techniques for 
three-dimensional imaging of target elements that overcome the problems of 
conventional three-dimensional imaging systems, including known tomosynthesis 
systems. The present invention enables the use of tomosynthesis to efficiently provide 
accurate three-dimensional imaging of a target element in a shorter time than has 
previously been possible. The invention includes an algorithm that is efficiently tuned 
by computationally appropriate segmentation of the target volume and directs the 

reconstruction of volume segments to a multiplicity of CPUs within a multi-CPU 

i 

computer cluster so that the volume segments can be reconstructed in parallel. In a 
preferred embodiment, the volume segments are selected to be optimally overlapping in 
order to provide mutual boundary coverage of what would otherwise be under- 
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determined edge volume segments. Thus it increases the speed of the overall volume 
reconstruction while providing a quality of reconstruction that is substantially the same 
as for reconstructions that do not employ the segmentation algorithm of the invention. 

5 The invention is applied below to one preferred embodiment in which the system is used 
for tomosynthesis mammography; however, the invention will be useful in a variety of 
three-dimensional imaging situations. For example, the invention can be applied to a 
variety of patient imaging problems such as heart imaging, or imaging of the soft tissues 
or bones of the hand. The imaging system of the invention can be used for diagnoses (as 

10 is described below for tomosynthesis mammography) or it may be used for other 

applications such as three-dimensional modeling for the purpose of fitting an implant 
(whether orthopedic, such as a hip or knee implant, an artificial heart, or other type of 
implant) or for use in surgical navigation systems. 

15 1. Introduction to an Exemplary Tomosynthesis Mammography System 

Tomosynthesis mammography is a three-dimensional breast imaging technique. 
It involves acquiring projection images of a breast at a plurality of viewpoints, typically 
over an arc or linear path. Three-dimensional distribution of x-ray attenuation coefficient 
of the breast volume is reconstructed from these projections. An exemplary 

20 Tomosynthesis system 10 of the invention for breast imaging is illustrated in Figure 1A. 
Tomosynthesis system 10 includes an X-ray radiation source or tube 12 which is 
generally directed toward a detector, in this example, digital detector 14. The tissue to 
be imaged, in this example compressed breast tissue 16 compressed by compression 
paddle 18, is located between X-ray source 12 and digital detector 14 so that the 

25 attenuation of a signal from the X-ray source and passing through the tissue can be 
measured by the detector. 

X-ray source 12 is rotatable about axis of rotation 20 in order to image 
compressed breast 16 from a plurality of angles. For this particular system 10, X-ray 
source 12 is located 44.3 em from axis of rotation 20 while the axis of rotation is located 

30 21 .7 cm above digital detector 14. In this embodiment, eleven projections are acquired 
by moving X-ray source 12 over a 50° arc (-25° to +25°) above breast 16 in 5° angular 
steps. Breast 16 and detector 14 are stationary during the image acquisition. While in 
this embodiment the X-ray source 12 is rotatable within a plane to different angles, it 
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should be understood that the source could be positioned so as to project images from a 
variety of viewpoints within the scope of the invention. 

Certain characteristics of this exemplary embodiment of a tomosynthesis system 
5 useful with the invention are described below: 

Spatial resolution and contrast resolution: Tomosynthesis system 10 can use an 
amorphous-Silicon-based flat panel detector 14 on which a Csl crystal phosphor is 
grown epitaxially read out as 2304X1800 pixels (lOOfim pixel pitch) via a TFT array. 
10 This particular detector has a linear response over exposure levels up to 4000mR and 12 
bits of working dynamic range. The reconstructed volume obtainable can be presented 
as a series of adjacent planes, with each plane having about the same resolution as the 
detector (lOOum), but with a depth resolution on the order of a millimeter. 

15 Dose: The target/filter combination is Rh/Rh and the accelerating potential is 25~33kVp 
to image breasts with 3~8cm range of thickness. The total x-ray dose for acquiring 1 1 
projections is approximately 1.5 times of that used for one film-screen mammogram. 
Each projection is a low dose breast image (approximately 1/1 1 of the does per 
projection). 

;20 

Patient motion: Patient motion is reduced by fast image acquisition. Using cone-beam x- 
ray geometry and an area detector, a projection of the whole breast can be recorded with 
one x-ray exposure at each angle. For each projection, the exposure time is 0.1~0.2s and 
detector readout time is about 0.3s. Rotation to the next angle is performed during the 
25 detector readout. The total image acquisition time for 1 1 projections is about 7sec. 
Breast compression also helps to reduce patient motion. 

Image acquisition geometry: The design of Tomosynthesis system 10 is based on the 
conventional mammography system. The MLO views have been used in most cases 
30 since it provides the most complete coverage of the whole breast. 
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Digital tomosynthesis mammography is particularly well suited for removing the 
tissue-overlap, and therefore reducing false-positive and false-negative diagnoses of 
breast cancer. The geometry of the tomosynthesis mammography system 10 is further 
5 illustrated in Figures IB and 1C which illustrate two orthogonal views of the 

tomosynthesis system geometry with the Figure IB view being along the patient's chest 
wall (with the X-ray source 12 traveling into and out of the page) and the Figure 1C 
view being in a direction toward the patient's chest wall. The system 10 acquires 1 1 
projection images of the breast and then reconstructs an estimate of the volumetric 

10 distribution of attenuation coefficients likely to have resulted in the measured 
projections.. Each projection image consists of 2304x1800 100 |um pixels. The 
thickness of compressed breasts ranges from <25 mm to >80mm, typically within this 
range. The chest-to-nipple distance ranges from <50mm to >160mm,but typically lies 
within this range. The reconstructed volume distribution consists can be represented by 

15 slices spaced 1 mm apart with lOOx 100 pm in-plane voxels. Therefore, the size of 
reconstruction image is (30-80) x 2304 x (500-1600). Using a PC with a 2.4GHz 
CPU, the reconstruction of a breast volume takes up to 5 hours. This is far to long for 
clinical use. One goal of this invention is to reduce the reconstruction time to a few 
minutes. This will not only permit clinical care, but will allow "real-time" needle 

20 placement to guide biopsies, etc. 

2. Volume Segmentation 

Projection images can be divided into a series of stripe-like segments, parallel or 
curvilinear oriented along the chest wall, from the chest wall to the nipple. Each 

25 projection segment covers part of the whole breast volume. With appropriate 

segmentation, projection segments will continuously inform the whole breast volume. 
The complete or partial volume can be composed by careful assembly and blending 
results of segment reconstructions. Increasing the number of segments results in smaller 
segmented volume coverage and faster reconstruction. 

30 Segmentation can be particularly useful in a parallel computing approach to 

image reconstruction. A parallel computing reconstruction can be developed and 
implemented on a computer cluster with, for example, 32-64 processors. In parallel 
computing reconstruction, the computation task is divided into independent smaller 
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tasks, each of which reconstructs a part of the whole breast accomplished by one 
processor of the computer cluster. The volume image of the whole breast can be 
retrieved from the results of the small tasks. In the Tomosynthesis methods and systems 
5 of the invention, reconstruction of segments can be divided among a number of 

processors so that the segment reconstructions can be performed in a simultaneous, 
parallel fashion. 

One method of segmentation, non-overlapping segmentation, can be illustrated 
by reference to Figure 2, which provides a virtual view that is from the same perspective 

10 as the view of Figure IB. The vertical axis is the breast thickness dimension (Z axis, 
typically between about 3 and 8 cm in thickness as marked on the Figure). The 
horizontal axis is the chest-to-nipple direction (Y axis, typically about 1 to 17 cm, shown 
as columns 100-1700 for the illustrated detector). The motion of the X-ray tube is 
generally perpendicular to this Z-Y plane (the plane of the drawing sheet), with 

15 intersection between the Z-Y plane and the X-ray tube occurring at a tube rotation of 0° 
(Z=66cm at 0° and Z=61cm at ±25°). In this example, the projection image is divided 
into continuous, non-overlapping segments (the segments being represented by the 
columns graphically illustrated in Figure 2) where each segment contains 100 columns. 
The volume coverage by each projection segment is shown by connecting the two edges 

20 of a segment (the bottom portion of each segment column illustrated in Figure 2) to the 
X-ray source (the right-most line of each pair of lines connects to the source at 0° (one 
example marked as element 24 and the left-most line of each pair of lines connects to the 
source at 25° (one example marked as element 26)). There is a mismatch between the 
volume coverage at a source rotation of 0° and that of a source rotation of 25°. This 

25 mismatch is bigger for thicker breasts and for segments further away from the chest wall 
as shown in Table 1 below which shows the mismatch of volume coverage at source 
locations of 0° and 25°. In the worst case (8cm thick breast, 16cm from the chest wall), 
this mismatch is -16 pixels. The mismatch of volume coverage indicates that, with this 
non-overlapping segmentation method, a small portion of the reconstructed volume (at 

30 the edge of the volume) will suffer from "missing projection data" within the segmented 
projection set. 
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Table 1. Mismatch of Volume Coverage by 0° and 25° Sources (with Pixel as Unit) 

In the above example of non-overlapping segmentation, a small region of "bad 
5 volumes" occurs at the edge of the reconstructed volume. In Figure 3 A, the areas 
marked by cross-hatching 28, 30 represent "good volumes" and the gap 32 between 
them represents the "bad volume". Overlapping segmentation can be used to improve 
this situation. In this method - overlapping segmentation — a projection segment, 
illustrated in Figure 3B, overlaps with its neighbor segments so that the volume 

10 coverage also overlaps. In Figure 3B, projection areas 28, 30 marked by cross-hatching 
now overlap in overlap region 34. If the overlap of projection segments is big enough, 
"bad volumes" 32 will be located only in the overlapping regions 34. Therefore, a "bad 
volume" region 32 for one segment 28, 30 may overlap with a "good volume" region in 
a neighboring segment. Data describing the entire reconstructed volume can then be 

15 retrieved by using only "good volumes" and ignoring redundant "bad volumes." 
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3. Preliminary Results: 50-Row Segments and 50% Overlapping 

Projection Segmentation: 

5 An example of the system and method of the invention was implemented using 

segmented reconstruction with 50-row segments and a high 50% overlapping. 
Projections were divided into multiple segments (1 to N+l) as illustrated in Figure 4 
from the chest wall edge 40 to the nipple edge with the chest wall at the bottom of the 
Figure (the geometry of Figure 4 is thus rotated 90 degrees counter-clockwise with 

10 respect to Figures IB, 2, 3A and 3B). Each segment consists of 50 detector rows. The 
first 25 rows of segment N overlapped with the last 25 rows of segment N-l and the last 
25 rows of segment N overlapped with the first 25 rows of segment N+l. Accordingly, 
as shown in Figure 5, Ni, the index of the first detector 14 row in a segment N, will 
equal 25 *N; N2, the index of the last detector 14 row in the segment N, will equal 

15 25*N+49; and N c , the index of the center detector 14 row in a segment N, will equal 
25*N+25. The total number of segment was determined by the size of the breast from 
chest wall to nipple. If the projection consists of M rows, the number of segments will 
be 2x(M/50)-l . All the eleven projections were segmented in the same way. 

The eleven projection segments with the same distance to the chest wall side 

20 were grouped into a set and used to reconstruct a volume segment. Reconstructed 
volume segments from all projection segments were then merged to form the whole 
breast volume. The reconstruction volume segment 50 from a projection segment had a 
"slanted rectangular" shape as shown in Figure 5. The slope of the "slanted rectangular 
volume" was determined by the location of projection segments used to reconstruct this 

25 volume. The slope of the volume can be represented by angle a as shown in Figure 5, 
where a is the angle made by the detector plane and the line connecting the center of 
segment N and X-ray source 14 located at a rotation of 0°. 

A reconstructed volume segment had "bad voxels" close to its boundaries 
because the mismatch of volume coverage by between projection segments at different 

30 angles as shown in Table 1 above. However, the consecutive projection segments were 
overlapped by 50%, so the boundary of one segment was close to the center of a 
neighbor segment. The corresponding reconstructed volumes overlapped in a similar 
way and this 50% overlap was sufficient to place all "bad voxels" in overlapping 
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regions. Therefore, a location taken by a "bad voxel" in one volume segment was 
overlapped by a "good voxel" in a neighbor volume segment so the whole breast volume 
can be retrieved with "good voxels 5 '. 

5 

Reconstruction of Segmented Volume: 

A segment of reconstruction volume was stored in a 3-D image array 52, as 
shown in Figure 6, consisting of parallel image slices 54. Because the volume had a 
"slanted rectangular" shape (the line connecting the centers of image slices 54 points to 

10 the X-ray source), there was a shift between slice centers to form array 52. The 

implementation of the reconstruction algorithm was almost the same as that for the 
conventional reconstruction method, except that an extra operation (shift of slice center) 
was now taken in the calculation from volume index (i, j, k) to coordinate position (x, y, 
z) and vice versa. The value of the shift was determined by the slice position (the 

15 distance above the detector and the slope of the slanted volume) and the slope of the 
slanted volume varied for each volume segment (represented angle a in Figure 5). 

Merging of Reconstructed Volume Segments: 

When all the volume segments were reconstructed, they were merged to form the 

20 image of the whole breast volume. Figure 7 shows the merging of two neighbor volume 
segments 56, 58. The two reconstructed segments had about 50% overlapping similar to 
the overlapping in projection segments. Bad voxels were in regions close to volume 
boundaries. However, the "bad voxel region" in one segment 60, 62 was covered by the 
"good voxel region" in a neighbor segment due to overlapping. The whole breast 

25 volume was then retrieved using only "good voxel regions". 

4. Performance of the Invention 

A Tomosynthesis patient case was selected for performance test of a segmented 
reconstruction method of the invention. The thickness of the breast was 7 cm and the 
30 distance from the chest wall to nipple was —14.5 cm. This was one of the larger target 
volumes from among those that had been imaged with the system of Figure 1 . 
Reconstructions were done with both a conventional (unsegmented) algorithm and the 
segmented algorithm for 9 iterations. The reconstruction volume was 2304x1450x70 
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(height, width and thickness) with a total size of -450MB. The workstation used for 
reconstruction had a Pentium III 866MHz processor and 4GB RAM. 

For the conventional algorithm, it took -53 min/iteration and a total of 8 hours 
5 for 9 iterations. For the segmented reconstruction algorithm, each projection was 

divided into 57 (2x[1450/50]-l) segments with a 50-row width. The first 25 rows of 
each segment were overlapped with the second 25 rows of the neighbor segment before 
it; and the second 25 rows of this segment was overlapped with the first 25 rows of the 
neighbor segment after it. The projection segments from different projections were 

10 grouped according to the distance from the chest wall. Therefore, 57 sets of projection 
segments were generated, each containing 1 1 projection segments with the same 
distance to the chest wail. Each set of the projection segments was used to reconstruct a 
volume segment. Each volume segment had 2304x50x70 voxels and took -19 minutes 
for the total nine iterations (~2.1min/iteration). Currently, the mergence of 

15 reconstruction segments takes -5 minutes for this case. These reconstructions were done 
one by one using a single computer processor so the total reconstruction time is the 
summation of time used for all segments. When implementing the algorithm of the 
invention on a computer cluster, the total reconstruction time will be mainly determined 
by the time for the reconstruction of one segment. 

20 Reconstructed volumes were compared. Four sets of reconstructed volumes 

were presented as slices (Z=10, 30, 50 and 70mm) and are shown in Figures 8 through 
1 1 illustrating these slices respectively. The images 8A, 9A, 1 OA and 1 1 A are from the 
conventional reconstruction algorithm, the images 8B, 9B, 10B and 1 IB are from the 
segmented reconstruction algorithm and the images 8C, 9C, 10C and 1 1C are the 

25 difference between the conventional reconstruction and the segmented reconstruction. 
The value of most pixels of the difference images is less than 1% of the value of the 
corresponding reconstruction image pixel. The result demonstrates that images from the 
segmented reconstruction algorithm provide equivalent diagnostic information as those 
from the conventional reconstruction algorithm. 

30 

5. Conclusion 

The segmented reconstruction method works successfully on a single processor 
system. In the test with a large-size target volume, the segmented reconstruction took 
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2.1 min/iteration for reconstruction with a 50-row projection segment and -5 minutes to 
merge the reconstruction segments to form the whole target volume. The image quality 
of volumes formed from the segmented reconstruction method is equivalent to that from 
the conventional reconstruction method. If this method were implemented on a 
computer cluster with the same CPU (Pentium HI, 866MHz, 133MHz bus speed, 
1 33MHz RAM speed) with a sufficient number of CPUs, the time for the whole 
reconstruction will be the same as the time for reconstructing one segment due to 
simultaneous, parallel processing. The total time for this case will be 

T = (2. 1 min/iteration) x (iterations) + 5min merging time + data transfer time 

The first two terms take -25 minutes for this case on the current computer. 
Compared with 8 hours conventional reconstruction, segmented reconstruction reduces 
the time by a factor of 20 while providing diagnostic quality images. The typical breast 
size of a patient is about the half of this testing case. 

Currently, computer clusters with 3 times faster processors (above 2.2GHLz, 
266MHz RAM, 533MHz bus speed) are available. With such a system, the total time 
(reconstruction + merging + transferring) for a case is expected to be reduced to —5 
minutes for normal breast size and <10 minutes for very large size. 

Advantages of the invention thus include fast computation speed for clinical 
applications with no loss of image quality. In addition, these performance increases can 
be provided with no increased cost at installations that already include multi-CPU 
computer clusters. Additionally, specialty FPGA arrays can be fabricated to implement 
multiple CPUs on a single PC add-in card. One or a series of these cards may employ 
the same strategy outlined above. 

A person of ordinary skill in the art will appreciate further features and 
advantages of the invention based on the above-described embodiments. For example, 
specific features from the references incorporated by reference above may be 
incorporated into systems, methods or computer program products of the invention as 
well as features referred to in the claims below which may be implemented by means 
described herein and described in those documents. Accordingly, the invention is not to 
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be limited by what has been particularly shown and described, except as indicated by the 
appended claims. 

What is claimed is: 
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CLAIMS 

1 . A tomosynthesis method for creating a three-dimensional reconstruction of a 
5 target element volume comprising: 

acquiring radiation absorbance images of the target element volume through a 
limited plurality of angles; 

dividing the target element volume into a plurality of volume segments; 
applying a reconstruction algorithm to each segment to generate a three- 
10 dimensional reconstruction of each volume segment; and 

merging the three-dimensional reconstruction of each volume segment to create a 
three-dimensional reconstruction of the target volume. 

2. The method of claim 1, wherein the Volume segments overlap. 

15 

3. The method of claim 2, wherein the volume segments overlap by an amount 
sufficient to result in a three-dimensional reconstruction of the target volume that does 
not differ substantially in quality from an unsegmented reconstruction. 

20 4. The method of claim 3, wherein the value of a majority of the pixels in the three- 
dimensional reconstruction of the target volume differ by less than about 1% from pixels 
in the unsegmented reconstruction. 

5. The method of claim 2, wherein the volume segments overlap by between about 
25 0 and 50 per cent. 

6. The method of claim 1, wherein the volume segments comprise volume 
segments having a complex shape that is dependent upon the acquisition geometry and 
makes the reconstruction of each segment independent of the reconstruction of any other 

30 segment. 
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7. The method of claim 1 , wherein the volume segments comprise volume 
segments having a slanted rectangular shape so that the reconstruction of each segment 
is independent of the reconstruction of any other segment. 

5 

8. The method of claim 1, wherein the images are obtained using an image 
acquisition element having: 

a radiation source positionable at a plurality of positions with respect to the target 
element; and 

10 a radiation detector positioned so as to detect radiation emitted by the radiation 

source passing through the target element and determine a plurality of attenuation value 
for radiation passing through the target element to establish a radiation absorbance 
projection image of the target element for a particular radiation source position. 

15 9. The method of claim 8, wherein the radiation source is positionable at a plurality 
of angles in a first plane. 

1 0. The method of claim 9, wherein the volume segments comprise volume 
segments having a slanted rectangular shape and having a base that corresponds to a 

20 plurality of consecutive detector pixel rows that are parallel to and spaced apart from the 
first plane, each slanted rectangular volume segment extending from its base in a 
direction toward the radiation source. 

1 1 . The method of claim 9, wherein the method is deployed in mammography. 

25 

12. The method of claim 11, wherein the first plane is substantially parallel to a 
patient's chest wall. 

13. The method of claim 12, wherein the volume segments comprise volume 

30 segments having a slanted rectangular shape and having a base that corresponds to a 

plurality of consecutive detector pixel rows that are parallel to and space apart from the 
patient's chest wall, each slanted rectangular volume segment extending from its base in 
a direction toward the radiation source. 
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14. The method of claim 1, wherein the method is carried out using a plurality of 
processors with at least one segment reconstruction being carried out using a first 
processor of the plurality of processors and at least one segment reconstruction being 

5 carried out using a second processor of the plurality of processors. 

15. The method of claim 14, wherein each segment reconstruction is carried out 
using a different processor of the plurality of processors. 

10 16. A system for three-dimensional tomosynthesis imaging of a target element 
comprising: 

an image acquisition element for obtaining a plurality of images of the target 
element from a plurality of angles having: 

a radiation source positionable at a plurality of positions with respect to the target 

15 element; and 

a radiation detector positioned so as to detect radiation emitted by the radiation 
source passing through the target element and determine a plurality of attenuation value 
for radiation passing through the target element to establish a radiation absorbance 
projection image of the target element for a particular radiation source position; and 

20 a processor configured to apply a reconstruction algorithm to the radiation 

absorbance projection images of the target element obtained from a plurality of radiation 
source angles to generate a three-dimensional reconstruction of the target element 
wherein the processor is further configured to divide the target volume into a plurality of 
/ image reconstruction volume segments for separate image reconstruction of the volume 

25 segments and mergence of the reconstructed volume segments into a three-dimensional 
reconstruction of the target element. 

17. The system of claim 16, wherein the system comprises a plurality of processors 
with at least one segment reconstruction being carried out using a first processor of the 
30 plurality of processors and at least one segment reconstruction being carried out using a 
second processor of the plurality of processors. 
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1 8. The method of claim 1 7, wherein the system comprises a number of processors 
that is at least equal to the number of image reconstruction volume segments that the 
system divides the target volume into and each segment reconstruction is carried out 

5 using a different processor of the plurality of processors. 

19. The system of claim 1 6, wherein the volume segments overlap. 

20. The system of claim 19, wherein the volume segments overlap by an amount 
10 sufficient to result in a three-dimensional reconstruction of the target volume that does 

not differ substantially in quality from an unsegmented reconstruction. 

2 1 . The system of claim 20, wherein the value of a majority of the pixels in the 
three-dimensional reconstruction of the target volume differ by less than about 1% from 

1 5 pixels in the unsegmented reconstruction. 

22. The system of claim 19, wherein the volume segments overlap by between about 
0 and 50 per cent. 

20 23. The system of claim 16, wherein the volume segments comprise volume 

segments having a complex shape that is dependent upon the acquisition geometry and 
makes the reconstruction of each segment independent of the reconstruction of any other 
segment. 

25 24. The system of claim 16, wherein the radiation source is positionable at a plurality 
of angles in a first plane with respect to the target element. 

25. The system of claim 24, wherein the volume segments comprise volume 
segments having a slanted rectangular shape so that the reconstruction of each segment 
30 is independent of the reconstruction of any other segment. 
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26. The system of claim 24, wherein the volume segments comprise volume 
segments having a slanted rectangular shape and having a base that corresponds to a 
plurality of consecutive detector pixel rows that are parallel to and spaced apart from the 

5 first plane, each slanted rectangular volume segment extending from its base in a 
direction toward the radiation source. 

27. The system of claim 24, wherein the system is configured to be deployed in 
mammography. 

10 

28. The system of claim 27, wherein the first plane is substantially parallel to a 
patient's chest wall. 

29. The system of claim 28, wherein the volume segments comprise volume 

1 5 segments having a slanted rectangular shape and having a base that corresponds to a 

plurality of consecutive detector pixel rows that are parallel to and space apart from the 
patient's chest wall, each slanted rectangular volume segment extending from its base in 
a direction toward the radiation source. 

20 30. A computer program for three-dimensional tomosynthesis imaging of a target 
element volume from a plurality of radiation absorbance projection images obtained 
from different positions from an image acquisition element having a radiation source 
positionable at a plurality of positions with respect to the target element and a radiation 
detector positioned so as to detect radiation emitted by the radiation source passing 

25 through the target element volume and determine a plurality of attenuation values for 
radiation passing through the target element to establish a radiation absorbance 
projection image of the target element volume for a particular radiation position, the 
computer program code being embodied in a computer readable medium and comprising 
computer program code for: 

30 dividing the target element volume into a plurality of volume segments; 

applying a reconstruction algorithm to each segment to generate a three- 
dimensional reconstruction of each volume segment; and 
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merging the three-dimensional reconstruction of each volume segment to create a 
three-dimensional reconstruction of the target volume. 

5 31. The computer program of claim 30, wherein the volume segments overlap. 

32. The computer program of claim 3 1, wherein the volume segments overlap by an 
amount sufficient to result in a three-dimensional reconstruction of the target volume 
that does not differ substantially in quality from an unsegmented reconstruction. 

10 

33. The computer program of claim 32, wherein the value of a majority of the pixels 
in the three-dimensional reconstruction of the target volume differ by less than about 1% 
from pixels in the unsegmented reconstruction. 

15 34. The computer program of claim 31, wherein the volume segments overlap by 
between about 0 and 50 per cent. 

35. The computer program of claim 30, wherein the volume segments comprise 
volume segments having a complex shape that is dependent upon the acquisition 

20 geometry and makes the reconstruction of each segment independent of the 
reconstruction of any other segment. 

36. The computer program of claim 30, wherein the radiation source is positionable 
at a plurality of angles in a first plane. 

25 

37. The computer program of claim 36, wherein the volume segments comprise 
volume segments having a slanted rectangular shape so that the reconstruction of each 
segment is independent of the reconstruction of any other segment. 



30 



WO 2005/055803 



PCT/US2004/039913 



-22- 

38. The computer program of claim 28, wherein the program code is configured to 
use images obtained using an image acquisition element having: 

a radiation source positionable at a plurality of angles in a first plane with respect 
5 to the target element; and 

a radiation detector positioned so as to detect radiation emitted by the radiation 
source passing through the target element and determine a plurality of attenuation value 
for radiation passing through the target element to establish a radiation absorbance 
projection image of the target element for a particular radiation source angle. 

10 

39. The computer program of claim 38, wherein the volume segments comprise 
volume segments having a slanted rectangular shape and having a base that corresponds 
to a plurality of consecutive detector pixel rows that are parallel to and spaced apart 
from the first plane, each slanted rectangular volume segment extending from its base in 

15 a direction toward the radiation source. ; 

40. The computer program of claim 38, wherein the program code is configured to 
be deployed in mammography. 

20 41. The computer program of claim 40, wherein the first plane is substantially 
parallel to a patient's chest wall. 

42. The computer program of claim 41, wherein the volume segments comprise 
volume segments having a slanted rectangular shape and having a base that corresponds 

25 to a plurality of consecutive detector pixel rows that are parallel to and space apart from 
the patient's chest wall, each slanted rectangular volume segment extending from its 
base in a direction toward the radiation source. 

43. The computer program of claim 30, wherein the program code is configured to 
30 be carried out using a plurality of processors with at least one segment reconstruction 

being carried out using a first processor of the plurality of processors and at least one 
segment reconstruction being carried out using a second processor of the plurality of 
processors. 
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44. The computer program of claim 43, wherein each segment reconstruction is 
carried out using a different processor of the plurality of processors. 



5 
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